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Abstract 



We present results on the unique reconstruction of a semi-infinite Jacobi operator 
from the spectra of the operator with two different boundary conditions. This is 
the discrete analogue of the Borg-Marchenko theorem for Schrodinger operators 
on the half-line. Furthermore, we give necessary and sufficient conditions for two 
real sequences to be the spectra of a Jacobi operator with different boundary 
conditions. 
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1 Introduction 



In the Hilbert space hi^) let us single out the dense subset //in(N) of sequences which have 
a finite number of non-zero elements. Consider the operator J defined for every / = {fk} 
in lfin{N) by means of the recurrence relation 



oo 
k=l 



{Jf)k ■■= h-ifk-1 + Qkfk + hfk+1 keN\{i} 

{Jf)i :=gi/i + &i/2, 



:i.2) 



where, for every n G N, 6„ is positive, while g„ is real. J is symmetric, therefore closable, 
and in the sequel we shall consider the closure of J and denote it by the same letter. 

Notice that we have defined the Jacobi operator J in such a way that 
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is the matrix representation of J with respect to the canonical basis in /2(N) (we refer the 
reader to [2] for a discussion on matrix representation of unbounded symmetric operators). 

It is known that the symmetric operator J has deficiency indices (1, 1) or (0, 0) [U Chap. 
4, Sec. 1.2] and [23l Corollary 2.9]. In the case (1,1) we can always define a linear set 
D{g) C Dom( J*) parametrized by (? G M U {+00} such that 

J* \ Dig) =: Jig) 

is a self-adjoint extension of J. Moreover, for any self-adjoint extension (von Neumann 
extension) J of J, there exists a 5^ G M U {+0x3} such that 

J(9) = J, 

[25| Lemma 2.20]. We shall show later (see the Appendix) that g defines a boundary condition 
at infinity. 

To simplify the notation, even in the case of deficiency indices (0, 0), we shall use Jig) to 
denote the operator J = J*. Thus, throughout the paper Jig) stands either for a self-adjoint 
extension of the nonself-adjoint operator J, uniquely determined by g, or for the self-adjoint 
operator J. 

In what follows we shall consider the inverse spectral problem for the self-adjoint operator 
J{9)- 
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It turns out that if J 7^ J* (the case of indices (1, 1)), then for all (7 G M U {+c>o} the 
Jacobi operator J{g) has discrete spectrum with eigenvalues of multiplicity one, i.e., the 
spectrum consists of eigenvalues of multiplicity one that can accumulate only at ±00, [25| 
Lemma 2.19]. Throughout this work we shall always require that the spectrum of J{g), 
denoted a{J{g)), be discrete, which is not an empty assumption only for the case J{g) = J. 
Notice that the discreteness of a{J{g)) implies that J{g) has to be unbounded. 

For the Jacobi operators J{g) one can define boundary conditions at the origin in complete 
analogy to those of the half-line Sturm-Liouville operator (see the Appendix). Different 
boundary conditions at the origin define different self-adjoint operators Jh{g), h G M U 
{+cxd}. Jo{g) corresponds to the Dirichlet boundary condition, while the operator Joo{g) 
has Neumann boundary condition. If J{g) has discrete spectrum, the same is true for Jh{g), 
V/i G M U {+C)o} (for the case of h finite see Section 2 and for h = 00, Section 4). 

In this work we prove that a Jacobi operator J{g) with discrete spectrum is uniquely 
determined by a{Jh^{g)), o'{Jh2{9)), with hi, /i2 G M and hi 7^ /12, and either hi or /i2. If hi, 
respectively, /i2 is given, the reconstruction method also gives /12, respectively, hi. Saying 
that J{g) is determined means that we can recover the matrix (11.31) and the boundary 
condition g at infinity, in the case of deficiency indices (1,1). We will also establish (the 
precise statement is in Theorem 13. 4p that if two infinite real sequences {Xk}k and {/ifcjfc that 
can accumulate only at ±00 satisfy 

a) {Afejfe and {fik}k interlace, i. e., between two elements of a sequence there is one and only 
one element of the other. Thus, we assume below that Xk < fik < Afc+i- 

b) The series Ylik^l^k — Xk) converges, so 

J]](Aifc - Afc) =: A < cx). 

k 

By b) the product TT r- is convergent, so define 

, _r Afe — An 

kj^n 

— I Xji -1 r fif^ Xji 



n 



^ , _r Afe — A„ 



c) The sequence {r„}n is such that, for m = 0, 1, 2, ... , 

\ 2m 

converges. 

Tk 

k 



d) For a sequence of complex numbers {/^fcjfc, such that the series 

k 



> converges 

^ Tk 
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and 

5:^ = 0, m = 0,l,2,... 

it must hold true that /3a; = for all k. 

Then, for any real number /ii, there exists a unique Jacobi operator J, a unique > hi, 
and if J 7^ J*, a unique 5^ G MU {+00}, such that a^Jh^i^g)) = {\k}k and a{Jh^{g)) = {fik}k- 
Moreover, we show that if the sequences {Xk}k and {fik}k are the spectra of a Jacobi operator 
J{g) with two different boundary conditions hi < h2 (/i2 G then a), b), c), d) hold for 
A = /i2 - hi. 

Necessary and sufficient conditions for two sequences to be the spectra of a Jacobi op- 
erator J{g) with Dirichlet and Neumann boundary conditions are also given. Conditions b) 
and c) differ in this case (see Section Hj). 

Our necessary and sufficient conditions give a characterization of the spectral data for 
our two spectra inverse problem. 

Our proofs are constructive and they give a method for the unique reconstruction of the 
operator J, the boundary condition at infinity, g, and either hi or /i2- 

The two-spectra inverse problem for Jacobi matrices has also been studied in several 
papers [lOl [El [151 (21! • There are also results on this problem in [9j. We shall comment on 
these results in the following sections. 

The problem that we solve here is the discrete analogue of the two-spectra inverse problem 
for Sturm-Liouville operators on the half-line. The classical result is the celebrated Borg- 
Marchenko theorem ^ [20] • Let us briefly explain this result. Consider the self-adjoint 
Schrodinger operator, 

Bf = -f"{x)+Q{x)f{x), xeR+, (1.4) 

where Q{x) is real- valued and locally integrable on [0, 00), and the following boundary con- 
dition at zero is satisfied, 

cosa/(0) + sina/'(0) = 0, aG[0,7r). 

Moreover, the boundary condition at infinity, if any, is considered fixed. Suppose that the 
spectrum is discrete for one (and then for all) a, and denote by {Afc(a)}fcgN the corresponding 
eigenvalues. 

The Borg-Marchenko theorem asserts that the sets {Afc(ai)}fcgN and {Afc(a2)}fceN for 
some ai 7^ 02 uniquely determine ai,a2, and Q. Thus, the differential expression and the 
boundary conditions are determined by two spectra. Other results here are the necessary 
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and sufficient conditions for a pair of sequences to be the eigenvalues of a Sturm-Liouville 
equation with different boundary conditions found by Levitan and Gasymov in (19] . 

Other settings for two-spectra inverse problems can be found in [3l U, Hi]. A resonance 
inverse problem for Jacobi matrices is considered in [7j. Recent local Borg-Marchenko results 
for Schrodinger operators and Jacobi matrices [131 12Z] are also related to the problem we 
discuss here. 

Jacobi matrices appear in several fields of quantum mechanics and condensed matter 
physics (see for example [5]). 

The paper is organized as follows. In Section 2 we present some preliminary results that 
we need. In Section 3 we prove our results of uniqueness, reconstruction, and necessary 
and sufficient conditions (characterization) in the case where hi and /i2 are real numbers. 
In Section 4 we obtain similar results for the Dirichlet and Neumann boundary conditions. 
Finally, in the Appendix we briefly describe -for the reader's convenience- how the boundary 
conditions are interpreted when J is considered as a difference operator. 



2 Preliminaries 



Let us denote by 7 the second order symmetric difference expression (see (11. ip . (11.21) ) such 
that 7 : / = {fk}km ^ {(7/)fc}fcGN, by 

i-ff)k:=bk-ifk-i + qkfk + bkfk+i, A;GN\{1}, (2.1) 
(7/)i:=gi/i + fei/2. (2.2) 

Then, it is proven in Section 1.1, Chapter 4 of [T] and in Theorem 2.7 of [23] that 

Dom(J*) = {fe kiN) : 7/ G hiN)}, J* f = lf, / ^ Dom(J*). 

The solution of the difference equation, 

(7/)=C/, CeC, (2.3) 

is uniquely determined if one gives /i = 1. For the elements of this solution the following 
notation is standard [1, Chap. 1, Sec. 2.1] 

P„_i(C):=/„, neN, 

where the polynomial -Pfc(C) (of degree k) is referred to as the k-th orthogonal polynomial of 
the ffist kind associated with the matrix (II. 3p . 
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The sequence {Pk{C)}'kLo ^/m(N) but it may happen that 

oo 

Y,\Pk{Of <oo, (2.4) 

fc=0 

in which case ( is an eigenvalue of J* and f{() the corresponding eigenvector. Since the 
eigenspace is always one-dimensional, the eigenvalue of J* is of multiplicity one . Moreover, 
since the (von Neumann) self-adjoint extensions of J, J{g), are restrictions of J*, it follows 
that the point spectrum of J{g), g E'RU {-|-oo}, has multiplicity one. 

The polynomials of the second kind {Qk{C)}'kLo associated with the matrix (11. 3p are 
defined as the solutions of 

bk^ifk-i + Qkfk + hkfk+i = Cfk, /c G N \ {1} , 
under the assumption that /i = and /2 = b^^. Then 

Qn-i(C) :=/«., neN. 

Qk{C) is a polynomial of degree k — 1. 

By construction the Jacobi operator J is a closed symmetric operator. It is well known, 
P, Chap. 4, Sec. 1.2] and [23], Corollary 2.9], that this operator has either deficiency indices 
(1, 1) or (0, 0). In terms of the polynomials of the first kind, J has deficiency indices (1, 1) 
when 

oo 

^|Pfe(C)|'<oo, forCGC\M 

k=0 

(this holds for all C ^ C \ M if and only if it holds for one ( G C \ M), and deficiency 
indices (0,0) otherwise. Since J is closed, deficiency indices (0,0) mean that J = J*. The 
symmetric operator J with deficiency indices (1, 1) has always self-adjoint extensions, which 
are restrictions of J*. When studying the self-adjoint extensions of J in a more general 
context the self-adjoint restrictions of J* are called von Neumann self-adjoint extensions of 
J [21123]. All self-adjoint extensions considered in this paper are von Neumann. 

Let us now introduce a convenient way of parametrizing the self-adjoint extensions of J 
in the nonself-adjoint case. We first define the Wronskian associated with J for any pair of 
sequences ip = {{pk}'^^i and ip = {V'fcjfcli hi.^) as follows 

Wkiip,^^) := bkiifki'k+i - 'ipky^k+i) , keN. 

Now, consider the sequences v{g) = {vk{g)}'^=i such that Wk E N 

Vk{g) ■■= Pk-i{0) + gQk-M , geR (2.5) 
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and 



(2.6) 



All the self-adjoint extensions J{g) of the nonself- adjoint operator J are restrictions of 
J* to the set [25l Lemma 2.20] 

Dig) ■■= {/ = {fk}T=i e Dom(J*) : lim Wn{v{g)j) = O} = 
= {fel2{N):^fel2{n),\imW4v{g)j)=0}. 

n— >oo 

Different values of g imply different self-adjoint extensions. If J is self-adjoint, we define 
J{g) := J, for all G M U {+cxo}; otherwise J{g) is a self-adjoint extension of J uniquely 
determined by g. We have defined the domains D{g) in such a way that defines a boundary 
condition at infinity (see the Appendix). 

It is worth mentioning that if J 7^ J* then, for all (? G M U {+00}, J[g) has discrete 
spectrum. This follows from the fact that the resolvent of J{g) turns out to be a Hilbert- 
Schmidt operator [251 Lemma 2.19]. 

Let us now define the self-adjoint operator Jh{g) by 

Jh{g) ■■= J{g) - h{-,ei)ei, heR, 

where {cfc}^]^ is the canonical basis in h^N) and (■,■) denotes the inner product in this 
space. Clearly, Jo{g) = J{g). 

We define Joo{g) as follows. First consider the space /2((2, C)o)) of square summable 
sequences {fn}^=2 aiid the sequence v{g) = {vk{g)}'kLi given by (12. 5p and (12.61) with g fixed. 
Let us denote Joo{g) the operator in /2((2, 00)) such that 

Joo{9)f = 7/, 

where (7/)^ is considered for any k > 2 and /i = in the definition of (7/)2, with domain 
given by 

DomiJM) ■■= {/ e /2((2,oo)) : 7/ G l2{i2,oc)), \imWn{v{g), f) = O} . 
Clearly, the matrix 

h2 






62 ■■■\ 

gs 63 ■ ■ ■ 

h 64 

64 gs '■• 
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which is our original matrix (11.31) with the first column and row removed, is the matrix 
representation of Joo(fi') with respect to the canonical basis in /2((2, oo)). 

It follows easily from the definition of Jh{g) that if J{g) has discrete spectrum, the same 
is true for Jh{g) (/i G M U {+00}). Indeed, for /i G R this is a consequence of the invariance 
of the essential spectrum -that is empty in our case- under a compact perturbation [22]. We 
shall show in Section H] that it is also true that Joo{g) has discrete spectrum provided that 
a{J{g)) is discrete. 

For the self-adjoint operator Jh{g), we can introduce the right-continuous resolution of 
the identity Ej^(g){t), such that Jh{g) = J^tdEj^^g^it). Let us define the function p(t) as 
follows: 

p{t) := (Ej,^(,)(t)ei,ei), tGR. (2.8) 
Consider the function (see [M] and [251 Chap. 2, Sec. 2.1]) 

mhiC, g) ■■= {{Jhig) - Ciy'e,, ei) , cxMg)) . (2.9) 

mfi{(,g) is called the Weyl m-function of Jh{g)- We shall use below the simplified notation 
m^Cg) := mo{(,g). The functions p{t) and m/i(C,fi') are related by the Stieltjes transform 
(also called Borel transform): 

dp{t) 

It follows from the definition that the Weyl m-function is a Herglotz function, i.e.. 



rrihiCg) 



fcvs,y; ^ ImC>0. 

ImC 

Using the Neumann expansion for the resolvent (cf. [251 Chap. 6, Sec. 6.1]) 

(Ms) - c/)-' = - J: + '^^^(Ms) - (ir , 

fc=0 ^ ^ 

where C G C \ (j{J{g)), one can easily obtain the following asymptotic formula 



mh{(,g) = -^ ^ ^5 + 0(C ), (2.10) 



as C ^ 00 (ImC > e, e > 0). 

An important result in the theory of Jacobi operators is the fact that m{(, g) completely 
determines J{g) (the same is of course true for the pair mh{C,g) and Jh{g)). There are two 
ways for recovering the operator from the Weyl m-function. One way consists in obtaining 
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first p(t) from m{(,g) by means of the inverse Stieltjes transform (cf. [25^ Appendix B]), 
namely, 

rb+5 

p(6) — p(a) = limlim — / {lmm{x + ie, g)) dx . 

The function p is such that all the moments of the corresponding measure are finite [H 
123]. Hence, all the elements of the sequence {t'^j^Q are in L2(M.,dp) and one can apply, 
in this Hilbert space, the Gram-Schmidt procedure of orthonormalization to the sequence 
{t'^j^Q. One, thus, obtains a sequence of polynomials {Pk{t)}'kLo normahzed and orthogonal 
in L2(M, dp). These polynomials satisfy a three term recurrence equation [23] 

tPk-iit) = h-iPk-2{t) + qkPk-lit) + hkPkit) ken\{l} (2.11) 
tPo{t)=qiPo{t)+hPi{t), (2.12) 

where all the coefficients 6^ {k G N) turn out to be positive and {k G N) are real numbers. 
The system (12.111) and (12.121) defines a matrix which is the matrix representation of J. We 
shall refer to this procedure for recovering J as the method of orthogonal polynomials. The 
other method for determining J from m{(,g) was developed in [12] (see also [21]). It is based 
on the asymptotic behavior of m{(,g) and the Ricatti equation 



bW''\C,9) = qn-C- , neN, (2.13) 

where Tn^'^\(,g) is the Weyl m-function of the Jacobi operator associated with the matrix 
(II. 3p with the first n columns and n rows removed. 

After obtaining the matrix representation of J, one can easily obtain the boundary con- 
dition at infinity which defines the domain of J{g) in the nonself-adjoint case. Indeed, take 
an eigenvalue. A, of J{g), i.e., A is a pole of m{(,g). Since the corresponding eigenvector 
fW = {fkWjkLi is in Dom{J{g)), it must be that 

hm Wn{v{g),f{\))=0. 

This implies that either lim„^oo W^n ({<5fc-i(0)}^]^, /(A)) = 0, which means that g = +oo, 

lim„^ooW^n(m-i(0)}r=i,/(A)) 



or 
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lim„^ooW^n({gfc-i(0)}r=i,/(A)) 



If the spectrum of Jh{g) is discrete, say a{Jh{g)) = {Xk}k, the function p{t) defined by 
(12. 8p can be written as follows 

p(t) = E-' 

' ' ni. 
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where the coefficients {ak}k are called the normalizing constants and are given by 



«n = ^|Pit(A„)|' . (2.14) 



fc=0 



Thus, ^/a^ equals the I2 norm of the eigenvector /(A„) := {Pk{Ki)}kLo corresponding to A„. 
The eigenvector /(A„) is normalized in such a way that /i(A„) = 1. 

Clearly, 

l = (ei,ei)= [ dp = J2 — - (2-15) 

Jr ^ Oik 

The Weyl m-function in this case is given by 

^ afe(Afe - C) 

From this we have that 

1 ^ A,-C ,1 



Therefore, 



a-' = lim (A„ - C)m(C, g) = - Res m(C, g) • (2.17) 



Let us now introduce an appropriate way for enumerating sequences that we shall use. 
Consider a pair of infinite real sequences {Xk}k and {iik}k that have no finite accumulation 
points and that interlace, i.e., between two elements of one sequence there is one and only 
one element of the other. We use M, a subset of Z to be defined below, for enumerating the 
sequences as follows 

VkeM Xk<l^k< Afe+i, (2.18) 

where 

a) If inffe{Afc}fe = -00 and sup;;,{Afc}fe = 00, 

M ■=Z and we require < < Ai . (2-19) 

b) If < supfc{Afc}fc < 00, 

M := {A;}fer-oo ' (^max > 1) and we require < < Ai . (2.20) 

c) If supfe{Afe}fe < 0, 

M := {A;}L-oo ■ (2.21) 
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d) If mfk{fik}k > 0, 



M := {A:}r=, 



(2.22) 



e) If -oo < inffc{/ifc}fc < 0, 

M := {fcj^fc^.^ , (fcmin < -1) and we require < < Ai . (2.23) 

Notice that, by this convention for enumeration, the only elements of {Xk}keM and {fik}keM 
allowed to be zero are Aq or ^q. 



3 Rank one perturbations with finite coupling con- 
stants 

In this section we consider a pair of operators JhAg) and J/ialfi'); where hi, G M, that is, 
rank one perturbations of the Jacobi operator J{g) with finite coupling constants. 

3.1 Recovering the matrix from two spectra 

Let (7 G M U {+C)o} be fixed. Since Jh{g) is a rank one perturbation of J{g)i the domain 
of J{g) coincides with the domain of Jh{g) for all G M. Moreover, since the perturbation 
is analytic in /i, the multiplicity-one eigenvalues, Xkih), and the corresponding eigenvectors, 
are analytic functions of h |16] . 

LEMMA 3.1. Let {Xk{h)}k be the set of eigenvalues of Jh{g) (h G 
following holds 

-TT^kih) = — , 

ah ak{h) 

where ak{h) is the normalizing constant corresponding to Xkih). 

Proof. For the sake of simplifying the formulae, we write Jh and X{h) instead of Jh{g) and 
Xkih), respectively {k is fixed). Let us denote by f{h) the eigenvector of Jh corresponding 
to X{h). Take any 5 > 0, taking into account that Dom(J/j+5) = Dom(J/i) and that Jh is 
symmetric for any G M, we have that 

^X(h + 5)-X{h)){f{h + 5),f{h)) = 
{Jh+sf{h + 6), f{h)) - {f{h + 6), Jhfih)) = 
= {{Jhu -Jh + Jh)f{h + 5), f{h)) - {f{h + 5), Jhf{h)) = 
= {{Jh+s-Jh)f{h + 5)J{h)) = -5. 



For a fixed k the 
(3.1) 
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Therefore, 



X{h + 6)-X{h) 1 1 

iim = — lim — 



5-0 6 s^o{f{h + 6)J{h)) ak{h)- 

□ 

The cornerstone of our analysis below is the Weyl m-function. Let us establish the 
relation between mh{(^,g) and m{C.,g). Consider the second resolvent identity 



{Mg) - CI)-' - {J{g) - CI)-' = (Jig) - Ci)-\J{9) - M9)){M9) - Ci)-' , (3.2) 

where C e C \ {cr{J{g)) U a{Jh{g))}. Then, for heR, 

mniC, g) - m(C, g) = {{(Jhig) - CI)-' - (Ag) - a)"') ei, ei) 

= {{Jig) - CI)-\H; e,)e,){Mg) - CI)-'e,, e,) 
= {h{{Jh{g) - CI)-'eu e,){J{g) - CI)-'e^, e,) 

= hmh{C,g)m{C,g) ■ 

Hence, 

REMARK 3.2. If J{g) has discrete spectrum, then mlCg) is meromorphic and, by fl3.3p . 
so is mh{C,g)- The poles of mh{Cig) are the eigenvalues of Jh{g)- Since the poles of the 
denominator and numerator in (13.31) coincide, assuming that h 0, the poles of mh{Ci g) cire 
given by the zeros of 1 — hm{C,g) and the zeros of mh{C,g) by the zeros of m{C,g). Thus, 
Jhiig) and Jh2{g) have different eigenvalues, provided that hi ^ h2. 

THEOREM 3.3. Consider the Jacobi operator J{g) with discrete spectrum. The sequences 
{/ifcjfc = (y{JhAg)) and {\k}k = cr{Jh2{g)),hi 7^ h2, together with hi (respectively, h2) 
uniquely determine the operator J, h2, (respectively, hi) and, if J J*, the boundary con- 
dition g at infinity. 

Proof. Without loss of generality we can assume that hi < h2. Consider the Weyl m-function 
m{C,g) of the operator J{g). Let us define the function 

ra{C,g) = '^''f/\ , CeC\R. (3.4) 

Notice first that the zeros of m{C,g) are the eigenvalues of JhAg) while the poles of m(C,5') 
are the eigenvalues of Jh2{g)- This follows from Remark 13.21 and (13. 4p . Let us now show that 
m{C,g) is a Herglotz or an anti-Herglotz function. Indeed, since m{C,g) is Herglotz, then 

1 — him ^ ^ —1 
h2m 



n^(C, g) = ^-1^ = 1 + , -1 - • (3.5) 



h2-hi {h2-hi)m{(,g) 
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Therefore, m(C, g) is Herglotz or anti-Herglotz depending on the sign of /i2 — hi. Recall that 
if a function / is Herglotz, then, — j is also Herglotz. Since h2 — hi > 0, m(C, g) is a Herglotz 
function. 

Thus, the zeros {fik}k of m{(,g) and its poles {Xk}k interlace. Let us use the convention 
f i2.l8p - fi2.23p for enumerating the zeros and poles of m{(,g). By this convention, if the 
sequence {\k}k (or {/^fc}fc) is bounded from below, the least of all zeros is greater than the 
least of all poles, while, if {Xk}k is bounded from above, the greatest of all poles is less than 
the greatest of all zeros. It is easy to verify, using for instance fl3.ip . that this is what we 
have for the zeros and poles of m(C, g) when J{g) is semi-bounded. 

According to flEl Chap. 7, Sec.l, Theorem 1], the meromorphic Herglotz function xn{(, g), 
with its zeros and poles enumerated as convened, can be written as follows 

"^(c-)-c^n'0-^)0-£)^ c>o, (3.6, 

where the prime in the infinite product means that it does not include the factor k = 0. 

From the asymptotic behavior of m{(,g), given by (12.101) . one easily obtains that, as 
C ^ 00 with ImC > e (e > 0), 

m(C, g) = l + {h- h,)C' + ih - h2)iqi - h2)C' + OiC') ■ (3.7) 

Therefore, 

lim m(C,^) = 1. 

ImC>e 

Then, using (13.61) . we have 

= hm rr' fl - ^Vl - ^) , e > 0. (3.8) 

Thus, m{(,g) is completely determined by the spectra a{Jh^{g)) and cr^Jh^ig))- Having 
found m{(,g), we can determine h2, respectively, hi, by means of (13. 7p . Hence, from (13.50 
one obtains m{C,,g) and, using the methods introduced in the preliminaries, J is uniquely 
determined. In the case when J ^ J* , we can also find the boundary condition g at infinity 
as indicated in Section 2. 

□ 

In [21] (see also [9]) it is proven that the discrete spectra of JhAg) and Jh^ig), together 
with hi and /12 uniquely determine J and the boundary condition g in the (1,1) case. Our 
result shows that it is not necessary to know both hi and /i2, one of them is enough. 
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It turns out that if one knows the spectra (j{Jh^{g)) and a{Jh,^{g)) together with gi, the 
first element of the matrix's main diagonal, it is possible to recover uniquely the matrix, the 
boundary conditions hi, h2 and the boundary condition at infinity, g, if any. Indeed, the 
term of order in the asymptotic expansion of m(C, g) 03. 7p determines hi — /i2. Since the 
coefficient of term is [hi — h2){qi — h2), if we know qi one finds /12, and then hi. 



3.2 Necessary and Sufficient conditions 

THEOREM 3.4. Given /ii G M and two infinite sequences of real numbers {Xk}k o-nd {^k}k 
without finite points of accumulation, there is a unique real /i2 > hi, a unique operator 
J{g), and if .J ^ J* also a unique G M U {+00}, such that, {iik}k = cr( Ails')) and 
{Xk}k = cr{'Jh2{9)) if and only if the following conditions are satisfied. 

a) {\k}k and {^k}k interlace and, if {Xk}k is bounded from below, mmk{fik}k > T^^T^k{Xk}k, 
while if {Xk}k is bounded from above, ma.Xk{Xk}k < ^SL,Xk{fik}k- So we use below the 
convention l^2.18i) - ^2.23\} for enumerating the sequences. 

b) The following series converges 

(/ife - Afc) = A < 00 . 

keM 



By condition [Jp the product TT ^ is convergent, so define 

Xk — An 

k^n 

:= ^^^^^ n ' VnGM. (3.9) 



A Xk — Ar 

keM 

kj^n 



c) The sequence {r^jneM is such that, for m = 0, 1, 2, . . . , the series 

\ 2m 

converges. 

Tk 



d) If a sequence of complex numbers {(3k}kGM is such that the series 

Ell^kf 
converges 

Tk 



and, for m = 0,1,2, ... , 

Tk 

k&M 

then l3k = for all keM. 



14 



Proof. We first prove that if {\k}k and {fik}k are tlie spectra of Jhiig) and Jh^{g), with 
/i2 > hi, then (H) , E]) , El) , and[^ hold true. The condition [aj) follows directly from the proof 
of the previous theorem. To prove that [6]) holds, observe that (13. ip implies 

f^'^ dh 



Consider a sequence {Mn}'^^i of subsets of M , such that M„ C M„+i and U„Af„ = M, then, 
using fl2.15p . we have 

The sequence {sn}'^=i is then convergent and clearly 

E(yUfc - Afe) = lim Sn = h2-hi. 
n— >oo 

Thus, A = h2 — hi. 

The convergence of the series in [6]) allows us to write (13.81) as follows 



c ^ = n':^ lim n' ^ 

l^k C~*oo J- J. At. - C 

Now, using again E]), it easily follows that for any e > 
Thus, C = TT /ifc/Afc and by (13. 6p . 



e > 0. 



m(C,^?)=nT^- (3-10) 



Let us now find formulae for the normalizing constants in terms of the sets of eigenvalues 
for different boundary conditions. By (12.171) . 

«n^(^2,^) = lim (A„ - Qmh^iCg) . 
Using the second resolvent identity, as we did to obtain (13. 3p . we have that 



1- {hi- h2)mh^{(,g) ' 
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Therefore, 



m(C, g) = ""'^i^'^l = l-{h,- h,)m,, , C e C \ M . (3.1i; 



Then, the normahzing constants are given by 

a-\h2,g) = lim (A„, - O^^^^l^l^ = j-^}im (A„ - (MCg) . 

C^A„ rl2 — ill ri2 — ''■I C^-^n 



Now, 



\--f(. / ' ' ' Xli 1 J / ^/ I I ^ u 

C-^-^n l^—rXn Ah — C 



hm (A„, - C)m(C, ^) = lim (A„ - C) JJ 



Hence, 



feGM " 

-t^n^- (3.13) 

^ ^ keM 

Notice that, since A = /i2 — /ii, it follows from (13.131) that r„ = a„ for all n G M. Hence 
the spectral function p of the self-adjoint extension Jh2{g) is given by the expression p(t) = 
'^x^:<t '^k^- ThusEj) follows from the fact that all the moments of p are finite [IlEB]- Similarly, 
stems from the density of polynomials in L2(M, dp), which takes place since p is iV-extremal 
P, [231 Proposition 4.15]. 

We now prove that conditions [oj) , Ej), [cj), and[^ are sufficient. Let {\k}k and {pk}k be 
sequences as indj) and [6]). Then, 

0<nf^<-- (3.14) 

k^n 

The convergence of this product allows us to define the sequence of numbers {Tn\n&M- Ob- 
serve that for all G M, r„ > 0. Indeed, A > and fl2A8D - ( 12:231 ) yield /i„ - A„ > for all 
n G M. Thus, taking into account (13.141) . we obtain 

r„>0, MneM. (3.15) 

Let us now define the function 

/^W-=E7' (3.16) 

Since (13.151) holds, p is a monotone non-decreasing function and has an infinite number of 
points of growth. Notice also that p is right continuous. Now, we want to show that for the 
measure corresponding to p all the moments are finite and 

dp{t) = l. (3.17) 
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The fact that the moments are finite follows directly from conditional). Indeed, 

r \m 

/ t-dp{t) = E ^ • 

We show next that (13.1 7p holds true. Given the sequences {\k}k and {iJik}k satisfying Ej) and 
[6]) we can define the function 

m(C):=nT-7- (3-18) 

Cell " ^ 

Taking into account (13.90 . one obtains that 

Res(m(C)-l) = --. 

C— An In 

In view of [6]), we easily find that 



lim (m(C) - 1) = hm J] " 1 

ImC>€ ImC>efceM ^ 



hm n fi+^^^Vi=o 



(3.19) 



f — >oo 

ImC>efc6A/ 



Thus, on the basis of Cebotarev's theorem on the representation of meromorphic Herglotz 
functions [iHl Chap. VII, Sec.l Theorem 2], one obtains 

A 



m(C) - 1 = E 7T-37wr • (3.20) 



fceM 



We now define the function m[Q := ^ . Then, (13.201) yields 



We next show that 



Indeed, 



lim C^(C) = ~1 

ImC>e 



A A J-1 Afc-C 



1 J 1 fl^k-C 

— exp 



.keM 



Xk-C 

P-k — Afc 



ie.p{Ei:(-ir'(t^)'' 

LfcGA/p=l ^ " ^ ^ , 
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Then, 



ImC>e 



lim C"^(C) 




Also, from (]3.2ip one has 



lim C"^(C) 



1 



Im f >e 



Therefore, 




Having found a function p with infinitely many growing points and such that fl3.17p is satisfied 
and all the moments exist, one can obtain, applying the method of orthogonal polynomials 
(see Section 2), a tridiagonal semi- infinite matrix. Let us denote by J the operator whose 
matrix representation is the obtained matrix. By what has been explained before, this 
operator is closed and symmetric. Now, define h2 ■= A + hi and J := J + /i2(-, ei) ei. 

If J = J*, we know that p{t) = jci, ei), where Ej{t) is the spectral decomposition 
of the self-adjoint Jacobi operator J. Then, obviously, J = Jh2- 

li J ^ J*, the Stieltjes transform of p is the Weyl m-function, we denote it by 'w{(), of 
some self-adjoint extension of J that we denote by J. This is true because of the density 
of polynomials in L2(M, dp). Indeed, means that the polynomials are dense in L2(M, dp). 
Thus, w{() lies on the Weyl circle, and then, it is the Weyl m-function of some self-adjoint 
extension of J pQ, [231 Proposition 4.15]. Therefore, J + h2{-, ei) ei is a self-adjoint extension 
of J and hence, J + /i2(-,ei)ei = J{g) for some unique (7 G M U {00}. Furthermore, we 
obviously have that, J = J/12 (5') ""^(C) = ''^h'zi.CiQ)- We uniquely reconstruct m{(,g) 
from mh^{C,g) using (13. 3p and then, we uniquely reconstruct g as explained in Section 2. 

Notice that we have 



It remains to show that a{Jh^{g)) = {Xk}k and a{.Jh^{g)) = {pk}k- To this end consider 
the function m{(,g) for the pair and J^-^: 




m(C,^) 



mh2{C,g) 
mh.iCg) ' 



C G C\M. 
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Let the sequence {'jk}k denote the spectrum of J/i^. Then, arguing as in the proof of (13.101) 
we obtain that 

7fc - C 



m(C,^7)= n 



Since we have aheady proven that [oj) and E]) are necessary conditions, we have that 

J^(7fc - Afe) = A < oo. 
fceM 

Then, as in the proof of (I3.19p . it follows that 

hm (m(C)-l) = 0. 

ImC>e 

Hence by Cebotarev's theorem [18', Chap. VII, Sec.l Theorem 2], 

h2 - hi 



where we compute the residues of m(C) as in fl3.12l) . Thus, since ak{h2,g) = t^, VA; G M, 

A 

(Afc - On 



fceAf 

But {Xk}k and {fik}k are the poles and zeros of xn{(,g) and then, the eigenvalues of Ajlfl') 
and Jhi{g), respectively. 

□ 

REMARK 3.5. We draw the reader's attention to the fact that the matrix associated with 
the function p, constructed in the proof of the previous theorem, may have deficiency indices 
(1,1) [Il[23l[25]. 

If we drop the condition of the density of polynomials in L2(ffi, dp) and our reconstruction 
method yields a nonself-adjoint operator J, then the sequences {Xk}k and {pk}k correspond 
to the spectra of some generalized self-adjoint extensions of and J/i^, respectively (see 
|23]). The generalized extensions of symmetric operators, which are not von Neumann 
extensions, were first introduced by Naimark (see Appendix I in [2J on Naimark's theory). 

In [15] the case of Jacobi operators bounded from below is considered. A uniqueness 
result is proven, and some sufficient conditions for a pair of sequences to be the spectra of a 
Jacobi operator with different boundary conditions are given. 
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4 Dirichlet-Neumann conditions 

4.1 Recovering the matrix from two spectra 

In this section we shall consider the pair of Jacobi operators Jo{g) = J{g) and Joo{g)- Here, 
as before, we keep the convention of writing J{g) even if J = J*. The matrix representation 
of Joo{g) corresponds to the matrix representation of J{g) with the first column and row 
removed. From the Ricatti equation (12.131) . taking into account that m^'^\C,^g) = m{(^,g) 
and m^^\C,g) = moo(C)5')) we have 



As before, we assume that the spectrum of J = J{g) is discrete. 

If m{(, g) is a meromorphic function, then, by (14.11) . moo(C; g) is also meromorphic and the 
spectrum of Joo{g) is discrete. The poles of m(^, g) are the eigenvalues of J{g), while the zeros 
of m(C, g) are the eigenvalues of Joo{g)- Since m{(, g) is always a Herglotz function, under our 
assumption on the discreteness of a{J{g)), m{(, g) is a meromorphic Herglotz function. This 
implies that a{J{g)) and (y{Joo{g)) are interlaced, that is, between two successive eigenvalues 
of one operator there is exactly one eigenvalue of the other. 

Let the sequence {Xk}k denote the eigenvalues of J{g) (the poles of m{(, g)). Furthermore, 
{/^fc}fc will stand for the eigenvalues of Joo{g) (the zeros of m{(,g)). It is worth remarking 
that, in contrast to the case of boundary conditions being rank one perturbations with finite 
coupling constant, here our convention for enumerating the sequences {Xk}k and {fik}k does 
not work in the case when J{g) is semi-bounded from above. Indeed, it follows from the 
mini-max principle pi] that if J{g) is bounded from below, the smallest of all poles is less 
than the smallest of all zeros of m{(,g), and if J{g) is bounded from above, the min-max 
principle applied to —J{g) implies that the greatest of all zeros is less than the greatest of 
all poles of m{(, g). 

So let us consider first the case when J{g) is not semi-bounded or semi-bounded from 
below and enumerate the sequences {\k}k and {/ifcjfc by (I2.18p . (I2.19p . (I2.22p . and (I2.23p . 
Then, by the same theorem we used to obtain (13. 6p |18j . m{(,g) can be written as follows 



where, as before, the prime in the infinite product means that it does not include the factor 




(4.1) 




(4.2) 



k = 0. 
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If J{g) is bounded from above, then we are still able to use (12.181) . (I2.20p and (12.211) for 
enumerating the zeros and poles of the meromorphic Herglotz function — . Thus, 



-^=C^n'fl-f)(l--)" , C>0. (4.3) 

Notice that, since we have enumerated zeros and poles of — by our convention, we have 
now 

WkeM, fik<Xk< /ifc+i , (4.4) 

and 

a) if < supfc{/ifc}fc < oo, 

M := {fcj^r-^^oo > (^max > 1) requiring A_i < < /n , (4.5) 

b) if supfc{/ifc}fc < 0, 

M:={k}t_^. (4.6) 

Here again Aq or fiQ are the only ones allowed to be zero. 
Equations (14. 2 p and (14. 3 P can be written in one formula 

where, if J{g) is not semi-bounded from above, K = C and {\k}k and {^k}k are enumerated 
by fl238D . (12191) . fl2:22|) . and fl2:23D . while = ~C~^ and {Aa,.}a,. and {/ifcjfc are enumerated 
by (I4.4l) - (l4.6p if J{g) is semi-bounded from above. 

We give now, for the reader's convenience, a simple proof of a theorem that was proven 
by Fu and Hochstadt [TU] for regular Jacobi operators (a regular Jacobi matrix is defined in 
[To]), and by Teschl [M] in the general case. 



THEOREM 4.1. (Fu and Hochstadt, Teschl) Consider the Jacobi operator J{g) with dis- 
crete spectrum. The sequences {Xk}k = '^{^{9)) (^iT'd {f^k}k = <^{Joo{g)) uniquely determine 
the operator J and, if J ^ J*, the boundary condition, g, at infinity. 

Proof. From (I2.10p we know that 

lim C^iCg) = -1, e > 0. 

ImC>e 
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Then, if J{g) is not semi-bounded from above, (14.2 1) yields 



where {Xk}k and {/ifcjfc are enumerated by flCTj) . fl239D . ([222]), and (12:231) . On the other 
hand, in the semi-bounded from above case (14.31) implies 

C-^ = Urn IY]' (l-^) (l-^) , e>0. (4.9) 

where {Xk}k and are enumerated by ( 14. 40 -( 1461) . Thus, in any case, one can find K, 

the constant in (14. 71) . from the sequences {Xk}k and {fik}k- Therefore, the spectra a{J{g)) 
and a{Joo{g)) uniquely determine m{(,g). Having found m{(,g) we can, using the methods 
introduced in Section [2l determine J and, in the case when J ^ J* , also find uniquely the 
boundary condition at infinity, g. 

□ 

REMARK 4.2. It turns out that, by (BSD and (jMl), K can be written as 

ir-^ = - hm Cn'fl--) fl-f) ' ^>0' (4-10) 

where the sequences {\k}k and {/ifcjfc have been enumerated by (12.181) . (I2.19p . (I2.22p . and 
( I2.23p . when J{g) is not semi-bounded from above and by (I4.4l) - (l4.6p . otherwise. 

In what follows the Weyl m- function will be written through (14.70 with K given by (I4.10p . 
From (14. 7p one can obtain straightforward formulae for the normalizing constants (I2.14p in 
terms of the sequences {Xk}k and {/ifcjfc. Indeed, when n 7^ 

hrn (A„ - C)m(C, g) = lim (A„ - 0^^ H = 

^ ^fcGM Afc 

=ir^(/.„-A„)^^n'^- 

Formulae (12.170 and (I4.10p then give 

r' / , \ / , \ -1 



«n^ = ^ 'f" ' / ' .-l' . ^^0. (4.11) 



22 



Analogously, 

-1 



^-1 ^ __ ^^keM\ 'J1Z_ . (4.12) 



ImC>e 



4.2 Necessary and Sufficient Conditions 

The following result establishes necessary and sufficient conditions for two given sequences 
of real numbers to be the spectra of J{g) and Joo(fi')- 

THEOREM 4.3. Given two infinite sequences of real numbers {Xk}k o-nd {^k}k without 
finite points of accumulation, there is a unique operator J{g), and if J ^ J* also a unique 
g G ]RU{+oo}, such that {\k}k = cr{J{g)) and {/Ufcjfc = o"(<^oo(fi')) if and only if the following 
conditions are satisfied. 

a) {\k}k o,nd {^k}k interlace and, if {\k}k is hounded from below, mmk{fik}k > m.mk{Xk}k, 
if {Xk}k is bounded from above, maxk{\k}k > niaxfc{/ifc}fc- So we use below the conven- 
tion Ii2.18\) . ^2.19\) . ^2.22i) . and Ii2.23\) for enumerating the sequences when J{g) is not 
semi-bounded from above, and ( [^.^[ j-( f7^ otherwise. 

By condition\W\j the product 



fceM 



converges uniformly on compact subsets of C (see the proof below and fT8\ Chap. 7, 
Sec.l]). 



b) The limit 



."-^n' (1-^)6-^)"' (-3) 



g6R keM ^ r / \ 



is finite and negative when the sequences {Xk}k and {fik}k are not bounded from above, 
and it is finite and positive otherwise. 

c) Let {r„}„gM be defined by 



- A„)^^ \\k^K, (\-^\(\-hL 

limg^ooi^TT' fl-f)f 



rn' = — w ' neM,n^O, 
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— 

The sequence {r„}„gAf is such that, for m = 0,1,2 ... , the series 

\ 2m 

converges. 



d) If a sequence of complex numbers {[3k}k&M, is such that the series 

|2 



1,3, 



k\ 



T~k 

k£M " 



converges 



and, for m = 0,1,2, ... , 



then /3k = for all k e M. 



J2^ = 0. 

keM 



Proof. We begin the proof by showing that the sequences a{J{g)) = {Xk}k and c^Jooig)) = 
{/ifcjfc satisfy dj), El), Ej), and[^. Since the Weyl m function is Herglotz, the eigenvalues of 
J{g) and Joo{g) interlace as indicated indj). To prove that E| holds, consider first the case 
when J{g) is not semi-bounded or only bounded from below, then (14. 81) yields E]). If J{g) is 
semi-bounded from above, (14 .yp implies [6]). 

On the basis of (14. lip and (I4.12p . r„ coincides with the normalizing constant a„ for all 
n G M. Hence the spectral function p of the self-adjoint extension J{g) is given by the 
expression p(t) = Yli\,,<t'^k^ ■ Thus [cj) follows from the fact that all the moments of p are 
finite [H [23]. Similarly, stems from the density of polynomials in L2(M, dp), which takes 
place since p is A^-extremal [Tj, [231 Proposition 4.15]. 

Let us now suppose that we are given two real sequences {\k}k and {pk}k that satisfy 
dj). It can be shown that 

k^n 

Indeed, the convergence of the infinite product follows from [aj) and is part of the Theorem 
1 in [18, Chap. 7, Sec.l] used to obtain (13. 6p . We give here, for the reader's convenience, 
some details. The product in (I4.14p converges if and only if 



An \ / . An \ . I . f 1 1 \ / . Ar ^ 



ky^n k^n 



< OO , 
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where prime means that the summand /c = is excluded. Thus, we have to prove that 



< oo . 



It will suffice to consider that in Ej) the sequences are ordered by fl2.18l) with M given by 
fl2A9l) . For any A; G N, (EUD implies 

Clearly, J2keN (l^ ~ convergent. Analogously, it can be proven that 



y(—-— 



< oo . 



fceN 

Having established the convergence of the the product in fl4.14p . its positivity follows easily. 

We have, therefore, a sequence of real numbers {r^j^gTy^- and let us now show that r„ > 0, 
G M. First notice that fl238D . fl2A9|l . fl2:22D . and ([223]), yield 

— (/in - ^n)^ — > (n 7^ 0) and /iq - Aq > . 
On the other hand (I4.4l) - fl4.6p imply 

— (/in - A„)^^^^^ — < (?i 7^ 0) and /io - Aq < . 

/in ^0 

From these last inequalities, taking into account fl4.14p and condition [&]) we obtain 

r„>0, VnGM. (4.15) 

Let us now define the function 

pW:=E-' (4.16) 

In view of f l4.15p . p is a monotone non-decreasing function and has an infinite number of 
points of growth. Now, we want to show that, for the measure corresponding to p, all the 
moments are finite and 

dp{t) = \. (4.17) 



•k 



The fact that the moments are finite follows directly from condition [cj). Indeed, 

AI 
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We show next that (14.171) is true. Given the sequences {\k}k and {fik}k satisfying dj) andE]), 
we can define the function 



m(C) : = 



1 



hnif 



keM 



1-ii 



Now, arguing as in the proof of (14.111) and (14.121) . we obtain 



Res m(C) = -r„ ^ 



On the other hand, 



n 



hm m(i^) = — hm 



fceA/ 



Afe 



keM 



1 



ii 



ii 



Thus, using again Cebotarev's theorem [18] we find that 

1 



fceAf 



7-fc(Afe - C) 



It follows from that 



lim i^m(i^) 

5— >oo 

Also from (14.191) one has 



Therefore, 



lim 

>oo . 



1-^ 



-1 



1-ii 



lim i^m(i^) 



fceA/ 



dp{t) . 



(4.18) 



(4.19) 



We have found a function p{t) with infinitely many growing points, such that all the moments 
exist for the corresponding measure and (14.171) holds. Therefore one can obtain, applying 
the method of orthogonal polynomials (see Section 2), a tridiagonal semi-infinite matrix. 
Let us denote by J the operator whose matrix representation is the obtained matrix. As 
was mentioned before, J is symmetric and closed. Now, if J = J*, we know that p{t) = 
{Ej(t)ei,ei), where Ej(t) is the spectral decomposition of the self-adjoint Jacobi operator 
J. If J ^ J*, then the Stieltjes transform of p{t) is the Weyl m-function m{(,g) of some 
self-adjoint extension of J with boundary conditions at infinity given by g, that is, 

dp{t) 



t-C 
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This last assertion is true because of the density of polynomials in L2{M.,dp), which follows 
from [3) • Hence p is A^-extremal [1]. This implies that m{(,g) lies on the Weyl circle, and 
then it is the Weyl m-function of some self-adjoint extension J{g) [1], |23j . 

It remains to show that a{J{g)) = {Xk}k and a{Joo{g)) = {pk}k- 

So we start from (I4.16p and find the Weyl m-function of J{g) using (14. 191) 



But {Xk}k and {pk}k are the poles and zeros of m and then the eigenvalues of J{g) and 
Joo{g), respectively. 

□ 

For Jacobi operators semi-bounded from below, necessary and sufficient conditions are 
given in p3] . Note that Remark 13.51 can also be made here. 

It is worth mentioning that, from (14.81) and (14. 9p . it follows that, when [6]) is seen as a 
necessary condition, one could write 

-1 



ImC>e 'sGM 

instead of (03]) . 



lim cTT'fi-^ 



Appendix: Boundary conditions for Jacobi operators 

The difference expression 7 defined by (12.1 p and (12. 2p can be written together in one equation 
with the help of some conditions. Indeed, consider the difference expression 7 given by 

{if)k = h-ifk-i + Qkfk + hfk+i , keN (60 = 1) • (A.i) 

Clearly, 7/ is equal to 7/ provided that 

/o = 0. (A.2) 

This requirement can be considered as a boundary condition for the difference equation 
flA.ip . Notice that, although /o is not an element of the sequence {/fc}^i, it can be used 
to introduce boundary conditions for flA.ip which turn out to be completely analogous to 



the boundary conditions at the origin for the Sturm-Liouville operator on the semi-axis. 
We shall refer to flA.2p as the Dirichlet boundary condition. Thus, J is the closure of the 
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operator which acts on sequences of /jj„(N) by (lA.ip with the Dirichlet boundary condition 
(lAj. 

Suppose that the deficiency indices of J are (1,1) and consider now the following solution 

of 

Vk{P) := Qk-M cos/? + Pk-M sin/3 , (3 e [0, vr) . 
Let us define the set 

{/ = {fk}k=i e km : 7/ G /2(N), hm W^{v{g)J) = O} . (A.3) 



Notice that D{g), defined by (12.71) . coincides with (]A.3|) as long as g = cot /3. As pointed out 



in Section 2, the domain of every self-adjoint extension of J is given by (1A.3P for some /3, 
and different /3's define different self-adjoint extensions [25]. Let us denote these self-adjoint 
extensions by J{g), as we did in Section 2, bearing in mind that g = cot /?. The condition 

lim Wn{v{g), /) = , / e Dom( J*) (A.4) 

n^oo 

determining the restriction of J* is considered to be a boundary condition at infinity. 

In analogy with the case of Sturm-Liouville operators one can define general boundary 
conditions at zero for the difference expression (lA.ip . To this end, consider the operator 
J{a,g) defined by the difference expression flA.ip with boundary condition at infinity flA.4l) 



if necessary, and boundary condition "at the origin" 

/i cosa + /o sina = , a G [0, vr) . (A. 5) 

Thus, if a G (0, vr), 

J{a,g) = J{g) - cota(-,ei)ei . 

Therefore J{a,g) = Jh{g), provided that h = cot a. 

When a = 0, from (lA.Sp . one has /i = and (lA.ip is used to define the action of the 
operator for k > 2. J{0,g) is said to be operator J{g) with Neumann boundary condition. 
For this case we have that J{0,g) is equal to Joo{g)- 

Acknowledgments We thank Rafael del Rio for a hint on the literature. 
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